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We analyze the thermodynamical properties of black holes in a modified theory of gravity, which 
was initially proposed to obtain correct dynamics of galaxies and galaxy clusters without dark 
matter. The thermodynamics of non-rotating and rotating black hole solutions resembles similar 
solutions in Einstein-Maxwell theory with the electric charge being replaced by a new mass depen¬ 
dent gravitational charge Q = ^/oGnM. This new mass dependent charge modifies the effective 
Newtonian constant from Gn to G = Gw(l + ce), and this in turn critically affects the thermo¬ 
dynamics of the black holes. We also investigate the thermodynamics of regular solutions, and 
explore the limiting case when no horizons forms. So, it is possible that the modified gravity can 
lead to the absence of black hole horizons in our universe. Finally, we analyze corrections to the 
thermodynamics of a non-rotating black hole and obtain the usual logarithmic correction term. 

PACS numbers: 


I. INTRODUCTION 

General relativity is one of the most successful and 
well-tested theories of gravitation to date. Nevertheless, 
there is still motivation to modify the framework at large 
scales. One of the main reasons is the discrepancy be¬ 
tween the observed dynamics of galaxies and clusters of 
galaxies and the amount of luminous matter these galax¬ 
ies and clusters contain mm- This is usually explained 
by postulating the existence of exotic dark matter. To 
date, no dark matter particle candidates have been de¬ 
tected in laboratory experiments or in satellite missions. 
An alternative resolution to the problem of the galaxy 
and galaxy cluster dynamics is a modification of the laws 
of gravitation on scales where Newtonian gravity or gen¬ 
eral relativity (GR) have not been extensively tested. 
One such framework called MOG (Modified Gravity) [4] 
has been able to explain the dynamics of galaxies and 
galaxy clusters without the need for dark matter in the 
present epoch of the universe [Ms]. 

In one version of the MOG formulation called the 
Scalar-Tensor-Vector Gravity (STVG) theory [4], the 
field content of general relativity has been increased to 
include scalar fields and a massive vector field. It has 
been used to describe the growth of structure, the matter 
power spectrum and the cosmic microwave background 
(GMB) acoustical power spectrum data in the early uni¬ 
verse [3 HO]. Solar system experiments are also in accor¬ 
dance with MOG m, and the dynamics of the Bullet 
Gluster has been explained without dark matter [12]. It 
may be noted that even though MOG fits the galaxy ro¬ 
tation curves like the phenomenological MOND ini, the 
advantage of using MOG is that it is constructed from 
a fully covariant modification of the Einstein-Hilbert ac¬ 
tion in GR. Thus, it is possible to apply MOG to cases 
where the gravitational field is strong, and the weak field 
approximation breaks down. It may be noted that there 
are various other ways to obtain the modified Newtonian 
potential from a covariant formalism. In fact, one such 


approach proposed by Bekenstein, called Tensor-Vector- 
Scalar (TeVeS) gravity produces MOND in the weak field 
approximation M- 

In this paper, we will analyze black hole thermodynam¬ 
ics in MOG. It may be noted that an entropy has to be 
associated with black holes in order to prevent the viola¬ 
tion of the second law of thermodynamics [mill]. Black 
holes also have a temperature and they evaporate due to 
Hawking radiation. However, it has been demonstrated 
that the evaporation of asymptotically de Sitter black 
holes (whose size is comparable to that of the cosmolog¬ 
ical horizon) differs significantly from the evaporation of 
asymptotically fiat black holes m- In fact, it has been 
observed that it is possible for quantum Schwarzschild-de 
Sitter black holes to anti-evaporate. The evaporation and 
anti-evaporation of black holes has also been analysed in 
/(T) gravity [18] and F{R) gravity [T9]- The black hole 
entropy gets generalized to Wald entropy [2Qll23] . The 
entropy formula can be used to compute the entropy of 
a black hole in any modified theory of gravity. The black 
hole entropy has been studied in Gauss-Bonnet gravity 
m [25]: Lovelock gravity [26l HZ] . The black hole 

thermodynamics has also been studied in massive grav¬ 
ity [28] 1^- 

The black holes have more entropy than any other ob¬ 
ject of the same volume [301EI], and this maximum en¬ 
tropy scales with the area of the black hole [32] . This ob¬ 
servation has led to the development of the holographic 
principle [33[E4]. However, as the black holes get smaller 
in size due to Hawking radiation, quantum fluctuations 
are expected to correct the standard relation between 
the area and the entropy of a black hole. This in turn 
can modify the holographic principle [351136] . Such cor¬ 
rections have been evaluated using several different ap¬ 
proaches. Non-perturbative quantum general relativity 
has been used to calculate the density of microstates for 
asymptotically fiat black holes m- Here, the density of 
states has been obtained by using a suitable conformal 
field theory. It has been demonstrated using the Gardy 



2 


formula that all black holes whose microscopic degrees 
of freedom are described by a conformal field theory will 
have logarithmic corrections [38l [39] . The exact partition 
function has been computed for BTZ black holes, and 
the logarithmic corrections have been obtained from an 
exact partition function [38|. Such correction have also 
been obtained using matter fields in backgrounds of a 
black hole [401442] . The logarithmic corrections are also 
produced from string theoretical effects |434[47| . Such 
corrections can also be generated from thermal fluctua¬ 
tions that occur in the thermodynamics of black holes 
[481149]. We will analyse the corrections due to thermal 
fluctuations for a MOG black hole, and observe that a 
MOG black hole also receives logarithmic corrections. 

In this paper, we will focus on static and rotating black 
hole solutions in MOG Uni ED, and will analyze their 
thermodynamics in the standard fashion. We show that 
the Hawking temperature and evaporation profile of the 
static MOG black hole is significantly modified over the 
standard Schwarzschild case. The rotating black hole 
is shown to admit remnant behavior. We also present 
corrections to our derived thermodynamic expressions, 
and find the usual logarithmic form. 


where 


- dycj)^. (7) 

The STVG field equations are given by |4]: 

G^v — + Qfii' = (8) 


where 


Q/x. = G(n0^^,-V^V,0), (9) 


where A is the cosmological con¬ 

stant and 0 = 1/G. 

The field equations for are given by 


= - J", 


( 10 ) 


and 


iiBya + — 0. (11) 

Moreover, we have the field equations 


II. MODIFIED GRAVITY ACTION AND FIELD 
EQUATIONS 

The action in the STVG formulation of MOG can be 
written as |4]: 


UG = K{x), 


njj, = L{x), 


where □ = V^Va and 


( 12 ) 

(13) 


S = Sg + S^ + Ss + Sm- (1) 


Here, So is the original Einstein-Hilbert gravity action, 
S(f) is the action of a massive vector field 0^, Ss is the 
action of scalar fields and Sm is the action for pressure¬ 
less matter: 


Sg 


Ss 


Sm 


^J^{R + 2A)V^gd‘x, 

I 


1C + V{ct>^) d^x, 


^ a0 f'C aG'CjjG Vq.^V;3^\ 

2^ \ J 


VciG) V^i^) 
G2 






( 2 ) 

(3) 

(4) 


■ j (PxAt^ + d‘^x + 


where R = g^^R^y, g = det{g^y)^ is the covariant 
derivative with respect to the metric g^y, the Vs are 
potential contributions to the action and we use units 
for which c = 1 and the metric signature is 
Moreover, 1C is the kinetic term for the field. It is 
possible to choose it to be the usual kinetic term for a 
vector field. 


K = 


(6) 



3 

87rG 


(l + 47r)V^GVaG 






and 


L(;r) = + ( 15 ) 


The covariant current density is defined to be 

Jf2 ^ (16) 

where is the energy-momentum tensor for matter 

and K = v^AGtv, = {G — Gn)/Gn is a scalar field, Gn 
is Newton’s constant, = dx^jds and s is the proper 
time along a particle trajectory. The perfect fluid energy- 
momentum tensor for matter is given by 

= {pM + Pm)u^u'' - pMg^\ (17) 


re Pm and pm are the density and pressure of matter, 
ectively. We get from (16) and (17) by using u^Uy = 


= SpMU^ • (18) 























3 


The gravitational source charge is given by 


Q = 


/ 


d^xJ^{x). 


(19) 


The values Q = and G = GAr(l + o^) origi¬ 

nate in the weak field approximation of the STVG field 
equations. The weak field approximation is based on a 
perturbation about the Minkowski metric r]^y\ 

9nv ~ 9iJiv T (^0) 


The test particle equation of motion is given by 
dPx^^ dx^ dx^ q dx^ 


( 21 ) 


where m and q = y/oG^m are the test particle mass and 
gravitational charge, respectively, and = 

1,2,3). Assuming that V^cj)^) is given by 


In our generalized gravitational theory electromag¬ 
netic waves (photons) and gravitational waves (gravi¬ 
tons) move with the speed of light [TT]. The range of 
values of the parameter a for weak gravitational fields 
is locally dependent on the size scale of the astrophysi- 
cal objects being investigated. A photon path is deter¬ 
mined by a null geodesic and the weak field solar system 
bending of light and the Shapiro time delay experimen¬ 
tal data are in agreement with GR for the parameter a 
having the bound a ^ 1 [4]. The perihelion advance of 
the planet Mercury’s orbit is in agreement with GR for 
(a ^ 1 [4]. The enhanced gravitational interaction expe¬ 
rienced by photons at the scale of galaxies and galactic 
clusters leads to an explanation of gravitational lensing 
data without dark matter for a > 1. 

In the following, we will use the matter-free MOG field 
equations [501 [5l]; 


= -SttGT^^^, (29) 


( 22 ) 

and djy(j)^ = 0, the weak field spherically symmetric 
static, source-free solution for fioir) is obtained from the 
equation (0 q = dfio/dr): 


= 0, (30) 

= 0. (31) 


00 + “00 “ /^^0o — 0. 
The solution is given by 


0o(r) = -Q 


exp(—/rr) 
r 


(23) 


(24) 


where the gravitational charge Q = ^/oGnM and M is 
the mass of the source particle. 

In the slow motion and weak field approximation, 
dr/ds ^ dr/dt and ^GMjr ^ 1, and for the radial ac¬ 
celeration of the test particle we get 


d^T GM 
dt‘^ 


qQ exp(-/ir) 
m 


(1 + /ir)- 


(25) 


For qQ/m = aGjsjM and G = GAr(l + <^), the modified 
Newtonian acceleration law for a point particle is given 

by 13; 


The energy-momentum tensor of matter in the 

gravitational field equations has been set equal to zero 
and we have 

= (32) 

We have set the cosmological constant A to zero and 
neglected the vector field mass ja. The best fits to 
galaxy and cluster dynamics yielded /i = 0.042 kpc~^ cor¬ 
responding to a mass = 2.6 x 10“^^ eV [5l|6]. This 
tiny mass can be neglected for the black holes solutions 
considered in the following. 

There is no reason to choose the usual kinetic term for 
the vector field 0^; it is possible to consider a non-linear 
kinetic term. Even though such a kinetic term will reduce 
to the usual kinetic term in the low energy limit, it can 
affect the behavior of black hole solutions. It is possible 
to consider general non-linear kinetic 1C terms [521458] : 


a{r) =-^—[I -ha — aexp(—/ir)(I + fir)]. (26) 

We observe that the acceleration of a particle is indepen¬ 
dent of its material content (weak equivalence principle). 
We can rewrite this modified acceleration equation as 


where the 
given by 

g{r) 


. , g{r)M 


(27) 


effective gravitational coupling strength is 


= GAr[I + a — aexp(—/ir)(I -h jir)]. 


(28) 


K = f{B), (33) 

where 

B = (34) 

A regular black hole soltion has been constructed by cou¬ 
pling the Einstein action to a non-linear electromagnetic 
field [52]. In fact, such a non-linear kinetic term has also 
been used in MOG for constructing regular black hole so¬ 
lutions [50| . We will explicitly consider such a non-linear 
kinetic term in MOG, without a singularity at radial co¬ 
ordinate r = 0. By using a Legendre transformation, we 
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can write /C = 2PHp — where dH = HpdP. Now for 
a regular black hole, we can write H as 


n 


P 


(l - 3^-2 
(l +v'-2Q2p)" 


3 ( V-2Q^P 

2Q^S yi+^-2 Q2p 


5/2 


(35) 


where s = QI2M = \JaG^ j2^ and P is a negative invari¬ 
ant quantity. The non-linear kinetic term corresponding 
to this can be written [52] : 


K = 


(l-8v/-2Q2p_6Q2pj 

(l + v/-2Q2p^^ 

3 (-2 Q2p)5/4 1^3 2^2 Q2p^ 


(l +v/-2Q2py^" 


• (36) 


By using these non-linear kinetic terms, we can derive a 
regular solution isni, without a singularity at radial co¬ 
ordinate r = 0. Such regular solutions in MOG based 
on nonlinear dynamics for the field can possess two 
horizons or no horizon, depending on the value of a EO]. 
For a < Q^crit = 0.673 a particular regular solution pos¬ 
sesses two horizons, while for a > acvit the solution does 
not possess horizons and is also regular at r = 0. The 
latter solution is called a “grey hole”. It is expected that 
Hawking’s information loss paradox [59l[6T] will not oc¬ 
cur for these gray holes. 

An electrically neutral black hole in MOG will appear 
similar to a GR static, spherically symmetric electrically 
charged Reissner-Nordstrom black hole [sniEi], but the 
MOG black hole solution depends only on the mass M 
and a. The rotating black hole solution in MOG will 
be similar to the Kerr-Newman solution in GR, but will 
depend on the mass M, the angular momentum J and not 
the electric charge. The Newtonian constant is modified 
by G = GAr(l + <^), where a is a free parameter. For weak 
gravitational fields the value of a is fixed by the fits of the 
modified Newtonian acceleration law to galaxy rotation 
curves and galaxy clusters: a = 8.89 laielE]. However, 
for the strong gravitational fields of MOG black holes, the 
weak field determination of a is no longer valid. For a 
given action for the vector field and the field strength 
the uniqueness theorems for the MOG solutions still 

hold. 

Real astrophysical bodies including black holes are not 
electrically charged. If they do have a small electric 
charge Qe, it will have a negligible effect on the curvature 
of spacetime [62]. For this reason the Reissner-Nordstrom 
black hole is not a physically interesting solution. On 
the other hand, the vector field in MOG is expected 
to be a physically significant field with positive gravi¬ 
tational charge Q = k,M = \/aG^M. Because there 


is no negative gravitational charge Q, MOG black holes 
are not gravitationally charge neutral. Furthermore, be¬ 
cause the gravitational charge in MOG is not indepen¬ 
dent of the mass of the black hole, it is not possible to 
construct static spherically symmetric extremal solutions 
where a remnant is left after Hawking radiation evapora¬ 
tion [30l[59]. In fact, in this paper it will be demonstrated 
that for static spherically symmetric MOG black holes, 
it is not possible to have a solution where the black hole 
temperature vanishes. 


III. THERMODYNAMICS OF STATIC 
NON-ROTATING BLACK HOLES IN MOG 


The first solution we will analyze is a static, spherically 
symmetric solution to the field equations. The metric for 
this solution can be written as uniEi]: 


ds^ 





de 



2GM 

-h 

r 



-1 


dr‘^ - r^dn^, 


(37) 


where the usual Newtonian constant Gw gets modified to 
G = Gw(l + and Q = y/oG^M. The metric exhibits 
two horizons: 


r± = GjsfM (l + (a ± \/l + <a) . (38) 

We observe that for a > 0, in contrast to the Reissner- 
Nordstrom solution when M < \Qe\ {Qe is the electric 
charge), the MOG black hole never has a naked singular¬ 
ity. Thus, as in the case of the Schwarzschild solution, 
the MOG black hole satisfies Penrose’s cosmic censorship 
postulate [63ll64) . 

The temperature of the black hole can be calculated 
via the surface gravity method evaluated at the outer 
horizon: 



1 dgoo 
2~^ 


(r = r+), 


(39) 


which gives 


2'kG]\[M 


1 

(1 + \/l + <^)(1 + <a + \/l + a) 


(40) 


This temperature can also be obtained by using the stan¬ 
dard expression 


r+ -r_ 
47rr^ 


(41) 


When a = 0, the usual Schwarzschild black hole tem¬ 
perature is obtained, T = I/^tiGnM. It is known that it 
is possible to construct an extremal solution for the elec¬ 
trically charged Reissner-Nordstrom black hole. When 





























5 


the temperature vanishes T = 0 a remnant exists af¬ 
ter the black hole has evaporated by Hawking radia¬ 
tion. Thus, the information loss paradox can be resolved 
for these extremal Reissner-Nordstrom black hole solu¬ 
tions. However, as we have explained, because electri¬ 
cally charged black holes are not expected to be phys¬ 
ically reasonable, this is not a viable resolution of the 
information loss paradox. However, as we will demon¬ 
strate, since an extremal solution does not exist for a 
static spherically symmetric MOG black hole as long as 
(a > 0, there is no solution for which T = 0 implying that 
there are no remnants remaining after Hawking evapora¬ 
tion. 

The temperature is plotted in Fig.[^ Solutions for a < 
0 are not considered, since it would result in a complex 
gravitational charge Q, an undefined temperature T as 
well as a negative gravitational constant. 



FIG. 1: Non-rotating MOG black hole Hawking temperature 
for a — 10,1,0.1 (from left to right). The gravitational con¬ 
stant has been set to Gw = 1- When o = 0, the temperature 
reduces to the Schwarzschild black hole temperature. Increas¬ 
ing values of the parameter a serve to suppress the radiative 
power for larger values of M, but with increasing power for 
M ^ 0. 


The radiative power of the black hole is 


— = AnrlaT\ (42) 

where a is the Stefan-Boltzmann constant. The lifetime 
is thus 


r 


f 

J Mbh 


dM 

dTTcrr^T^ 


Figure shows the ratio of the lifetime for various val¬ 
ues of a to that of a Schwarzschild black hole (a = 0). 
The dependence on a clearly influences the lifetime of 
such black holes. A value of a = 2 increases the value 
by roughly a factor of 20. A two-fold increase can be 
obtained with as small a value as a = 0.312. 



FIG. 2: MOG to Schwarzschild black hole lifetime ratios t/tq 
for various values of a. The gravitational constant and the 
Stefan-Boltzmann constant have been set to Gn — ct — 1. 


Although algebraically similar to the Reissner- 
Nordstrom black hole solution, the MOG metric is fun¬ 
damentally different in that the lack of electric charge 
implies that no zero-temperature remnants may be ob¬ 
tained (in contrast to the extremal Reissner-Nordstrom 
case). This is due to the fact that the gravitational charge 
Q is proportional to M. We note the “extremal limit” 
= r_ occurs only when a = —1, at which point the 
horizons vanish. 

The entropy of the black hole can be derived from the 
temperature by 


S'Y = 


Jo 


dM 


— TrGwAfgjj (\/l “h Ck. + 1 + 0^) (a/I T ol + l) . (44) 

Written in this form, it is apparent that it is related to 
the entropy of the Schwarzschild black hole as given by 
the area-entropy law. Defining A+ = 47rrl and using 
Eq. dM, we find 


5A = ^=7rGjvM|H(l + \/r+^)'. (45) 


47r^ 

3(7 


(1 -|- \/l T ot \/ir-l-~^)^GwAfgjj 


(43) 


Thus, the entropy defined in standard thermodynam¬ 
ics is different from the entropy calculated using the 
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Bekenstein-Hawking bound [MISO]. This can be taken 
as an indication that the Bekenstein-Hawking bound 
is modified in MOG. The MOG correction to the 
Bekenstein-Hawking bound can be written as 

AaS = St ~ Sa — ttG jsf + VlT-F^). (46) 

The MOG leads to an interesting modification of black 
hole thermodynamics. The entropy reduces to the famil¬ 
iar Schwarzschild black hole value S = when 

(a = 0. 

The local and global thermodynamic stability of a 
black hole is determined by both the heat capacity and 
the free energy, respectively IS3IM]. A calculation of the 
heat capacity for the MOG black hole yields: 

C =-27rGivM2(l + \/rTa)(l+a+\/rTa) < 0 , (47) 

which typically indicates the black hole is thermodynam¬ 
ically unstable for all values of M. The free energy F can 
be obtained from 


F = M -TS (48) 

Note the result will differ depending on the choice of the 
entropy expression. Ghoosing the thermodynamic defini¬ 
tion ( |HI[ ), one finds 


F=^ 

2 


while choosing the area entropy (45) yields 


(49) 


F = M[1-^ 


1 + \/l + (a) 
1 + (a H" ^/T^ 


a 


(50) 


using (HI). Stability of the thermal system occurs when 
F is minimized. For (49), this is at M = 0, and so 


the system is globally unstable. In the latter case (50), 
M = 0 is again the minimized state since F > 0 for 
Of > 0. The MOG black hole is thus both locally and 
globally thermodynamically unstable. 


IV. THERMODYNAMICS OF ROTATING 
BLACK HOLES IN MOG 

A solution for the rotating black holes in MOG has 
been constructed [50]. We will call this the Kerr-MOG 
solution, as it reduces to the usual Kerr solution in the 
limit Of = 0. The Kerr-MOG metric can be written as 
(a = J/M): 

ds^ = ^ (dt — a sin^ 0 dcjSj ^ ^ [(r^ + a^)d(j) — adt\ ^ 

— ^dr^—p^dO^^ (51) 


which describes black holes with horizon radii: 


T T — G 7 v(lTQ^)Af 


1 ± 


I 0 ^ a 

G%{lFa)‘^M‘^ IF a 


(52) 

Moreover, A = — 2GMr F o? F aGNGM‘^ and = 

+ oS cos^ 0. The ergosphere is located at 


Te = GAr(l+Q^)M 


1 + W1 — 


cos^ 0 


a 


G%{lFayM‘^ IF a 

(53) 

This solution is again algebraically identical to the Kerr- 
Newman metric but now with the gravitational charge 
Q = y/oG^M. The temperature for this Kerr-MOG 
black hole solution is determined by 


T = — 
2tt 


n = 


We obtain 


T = 


1 


— r_ 

2 {r\ F o?) 




2t^G jAFI 2^{oi -|- 1) T T 2 


(54) 


(55) 


where 


A- 


1 - 




\ (l + a)2 


(56) 


In contrast to the spherically symmetric static MOG 
black hole, the Kerr-MOG back holes possesses an ex¬ 
tremal configuration when , or 


1 - 


a 


G%M‘^{lFay IF a 


= 0 


ALxt — 


Gn \/l + 


a 


(57) 


The Kerr-MOG black hole free energy can be calcu¬ 
lated from (48), using the temperature ( [^ and entropy 
defined by (HI), and is displayed in Figure The free 


energy is minimized at the same point that defines the 
transition point for the heat capacity. This strengthens 
the assertion that this represents a critical point in the 
thermodynamic stability, and the positivity of the free 
energy indicates global instability so the black hole un¬ 
dergoes complete evaporation. 

It can be shown that the temperature (55) vanishes 


for a value of a, implying that this is a zero-temperature 
remnant. Gorresponding temperature profiles are plotted 
in Figure and the generic behavior of the heat capacity 
is shown in Figure The effect for increasing a > 0 is to 
shift the temperature curve to the left, resulting in small 
remnants. Furthermore, the maximum temperature of 
the Kerr-MOG black hole is correspondingly decreased. 
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FIG. 3: Kerr-MOG temperature profile for spin parameter 
a = 1 and a = 0,0.5,1. The remnant mass values for 
each curve are Mq = 1,0.816,0.707 for a = 0,0.5,1, re¬ 
spectively. These reach maximum temperatures of Tmax = 
0.023,0.019,0.016 for M^ax = 1.272,1.040,0.902, respec¬ 
tively. The curves do not appear to end at Th = 0 due to 
plotting resolution. 


V. REGULAR MOG BLACK HOLE 

In this section, we will analyze the thermodynamics of 
regular solutions which occur in the full non-linear MOG 
action. It may be noted that a regular solution for electri¬ 
cally charged objects in GR has been constructed using 
non-linear electrodynamics [52]. As the kinetic part of 
the action of the vector field in MOG can be taken to 
be non-linear, it is possible to construct a regular black 
hole solution in MOG [50]. This MOG regular solution 
can have two horizons, an outer horizon and an inner 
Cauchy horizon r_ for a < acrit = 0.673. However, for 
a > acrit there are no horizons and no naked singularity 
exists. The case without a horizon is called a “gray hole”, 
and the gravity of this gray hole for a sufficiently large 
mass M is expected to be so strong that it effectively 
can be treated as a black hole. The redshift observed by 
an asymptotic external observer will be large and there 
is small probability for any object to leave its gravita¬ 
tional field. However, as this probability is not zero, it 
will not suffer from the black hole information loss para¬ 
dox [59H6T]. So, it is interesting to note that the modifi¬ 
cation of gravity which was phenomenologically proposed 
to explain the correct dynamics of galaxies might end up 
predicting the absence of a horizons for objects in our 
universe, and hence end up solving the black hole infor¬ 
mation loss paradox for these astrophysical objects. 

The metric function for the regular solution can be 


FIG. 4: Kerr-MOG heat capacity for 0=1 and a = 1. The 
black hole is thermodynamically unstable when C < 0, and 
reaches a maximum temperature when C ^ ±oo (Mmax ~ 0.9 
in this case). For M < Mmax, the heat capacity is positive 
and the black hole evaporates to a locally stable cold remnant. 


written as [50] : 

dr^ 

ds^ = f{r)dt^ - (58) 

where 

_ 2GMr^ aGGNM\^ 

^ (r2 + aGNGM^f^^ ^ + aGGNM'^f' 

(59) 

It describes a Schwarzschild-MOG black hole for large r 
and the metric becomes asymptotically fiat in the limit 
r ^ oo. For small r, it describes a de Sitter (anti-de 
Sitter) core depending on the value of a with the effective 
cosmological constant: 

, 3 f 0^/2 _ 2 \ 

~ G%M^ \ay^{l + a)) ' 

The metric function /(r) approaches 1 in the limit r = 0. 
Horizons are obtained as the roots of /(r) = 0, which 
depend on the critical value a = acrit = 0.673. Figure 
shows the behavior of /(r) for a normalized mass M = 1 
and Gn = 1, in which the three phases of solutions are 
visible. 

A gray hole can approach a black hole with horizons 
as a limiting process, when a acrit- In general, the 
thermodynamic properties of the gray hole cannot be the 
same as a MOG black hole possessing horizons. It is 
expected that a gray hole’s entropy will be less than the 
entropy of a black hole of a similar size and mass, for 
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FIG. 5: Kerr-MOG Gibbs free energy for a — 1 and a — 1. 
The curve reaches a minimum near M ~ 0.9, consistent with 
the transition point of the heat capacity. 



FIG. 6: Regular MOG metric function for values a = 0.3 
(lower curve), a = acrit ~ 0.673 (middle curve), and a = 1 
(upper curve). The black hole mass is normalized to M = 1, 
as is the gravitational constant Gw = 1. The radial value for 
which /(rcrit) = 0 is rcrit ~ 1.68. Note that the value of acrit 
is independent of the mass. 



it is not a maximum entropy object. It is possible to 
study the thermodynamic features of the regular solution 
as a ^ Q^crit, when horizons form and an asymptotic 
observer will see that the redshift of light will approach 
infinity. 

If MOG is the correct theory describing gravity, then 
it is possible that only a regular solution without hori¬ 
zons forms during gravitational collapse. Because there 
are no horizons the gray hole will have negligible or no 
Hawking radiation and there is no information loss para¬ 
dox. However, we must caution that at present, we have 
no known physical law that demands that a > (acrit- 
In order to evaluate the temperature for the regular 
black hole with a < Ofcrit: we proceed by the surface 
gravity method. The derivative of the metric function is 
given by 

, 4(1+ a)Mr 6(l + (a)Mr^ 

/ (^) =-^ 37 ^- + -^ 57 ^- 

2a{la)M‘^r 4(a(l + (a)M^r^ 

where C = r^ + (a(l + a)M‘^ and the derivative is eval¬ 
uated at the outer horizon r = The corresponding 
real, positive roots of /(r) = 0 are given in Table Eval¬ 
uating the temperature T = /'(r = r+)/47r for the values 
therein, one observes that T ^ 1/M for all values of r. 
This means that in the presence of horizons, the tem¬ 
perature diverges as the black hole evaporates, despite 
the regular nature of the solution at r = 0. Thus, the 
divergent behavior of the temperature T in the singu¬ 
lar Schwarzschild black hole solution is retained for the 
regular black hole solution with an outer horizon. 


a T- r+ 

0.1 0.163M 2.067M 
0.2 0.135M 2.117M 
0.3 0.438M 2.148M 
0.4 0.676M 2.152M 
0.5 0.904M 2.120M 

TABLE I: Inner and outer horizons r± for the regular MOG 
black hole for various values of o > 0. 


VI. MOG BLACK HOLE ENTROPY 
CORRECTIONS 

It is expected that quantum fluctuations in the geome¬ 
try will lead to thermal fluctuations in the thermodynam¬ 
ics of black holes. The explicit form of such corrections 
can be calculated by analyzing the thermal fluctuations 
around a state of equilibrium [48l [49] . This can be done 
using the path integral formalism. In the path integral 
formalism, it is possible to calculate the amplitude for a 
field configuration to propagate to another field configu¬ 
ration. This can be done using the Euclidean quantum 
gravity formalism, where the temporal coordinate is ro¬ 
tated in the complex plane isTHn]. Thus, the partition 
function for spacetime can be written as 

Z = j [D]eM-SE), (62) 

where Se = —iS is the Euclidean action corresponding 
to the MOG action S. Thus, Se is obtained from S by 
a rotation of the time axis in the complex plane. 
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It may be noted that the partition function can 
be related to the statistical mechanical partition func¬ 
tion uaiTa]: 


poo 

= / dEp{E)exp{-/3E), 

Jo 


(63) 


where p is the inverse of the temperature. The density 
of states can be written as 


1 pPo^ 

^ / 

JBn-v 


/3o+ioo 


= 7 WW / d^exp[S'(/3)], 


where 


Po—ioo 


S = i3E-^\nZ. 


(64) 


(65) 


Usually this entropy is measured around the equilib¬ 
rium temperature /do, and all thermal fluctuations are 
neglected. This is done by making the identification 
T = /d“^, and Sq = S{/3)p=pQ = AJAGn, where as before 
A is the area of the black hole m- Identifying the tem¬ 
perature in (|^) with /do and using the entropy expres¬ 
sions (III) and (45) for S{/3o) and S'(/do), respectively, we 
can write the difference: 


A5(/3o) = 5(^o) - ^(/3o) 

= ttG ]\[T \/1 T^). 


( 66 ) 


Taking such thermal fluctuations into account, we can 
expand S{/3) around the equilibrium temperature /3o [48l 
|49] to obtain 




5 = 5o + -(/3-/3o)^(^^j 


(67) 


/ 3 =/ 3 o 


where we have neglected higher order corrections to the 
entropy and defined 5'o = S{P)\p=pQ. The density of 
states can now be written as 


p{E) = 


exp(5'o) 

27ri 


i 


PoB-ioo 


dp 


^Q—ioo 


1 


fd^S{p)\ 


xexpl-(/?-/?o)^^ 


)s- 


( 68 ) 


/ 3 =/ 3 o> 


After a change of variables, we get 
exp(S'o 


and it follows that 


5 = 5o - - In 


d^m'] 
dP^ A= 


/ 3 =/ 3 o 


- 1/2 


V dp^ T= 


/ 3 =/ 3 o 


(69) 


(70) 


The second derivative of entropy is actually a fluctuation 
squared of the energy. 


It is possible to simplify this expression by using the 
microscopic degrees of freedom calculated from a confor¬ 
mal field theory. The modular invariance of the confor¬ 
mal field’s partition function constraints take the form 
5'(/d) = U/d + b/3~^ [38], which can be generalized to the 
expression: S{/3) = a/3^ -h [48l [49], where all the 

constants satisfy, m^n^a^h > 0. This has an extremum 
at /do = (n5/ma)^/"^+’^ = T~^. 

Expanding the entropy around this extremum, we get 

S{P) = [(n/m)™/(’"+”) + (m/n)”/(’"+”ya”6™)i/(’"+”) 

+ i[(m + n)TO(”+2)/(™+”)n(’"-2)/(m+n)] 

x(„n+2/,m-2)l/(m+n)(/j_^^)2^ (71) 

By comparing Eq. ( |67| to Eq. , we obtain 

= (n/m)™/(’"+”) + (TO/n)”/(™+”)(a”6™)i/(™+”\ 

(72) 

and thus 

\ op J .a 


^ ^^n+2^m-2^1/(m+n) 


(73) 


We can now find the values of and simplify the ex¬ 
pression to obtain: 


V dp^ 


= ymnSoT^ 


(74) 


/ 3 =/ 3 o 


where 

3^mri 


^ ^^^(n+2)/(m+n)^(m-2)/(m+n) \ 
(n/m)W(7n+n) (m/n)V(7n+n) J 


r fi X 2/(m+n)' 

ktu) 


(75) 


As the factors ymn are independent of any black hole 
parameters, they can be absorbed by a suitable redefini¬ 
tion (see iiHiiin])- We will use the corrected form for the 
entropy, which neglecting higher order corrections can be 
written as 


S = So-]^\nSoT\ 

We can now write the the corrected entropy as 

Sa = 7rG'ivM2(l +vTT^)2 


(76) 


1 


1 


-1 In ■ _ . 

2 47rGAr (1 + a + Vl + a)"^ 


(77) 


This is the expression for the corrected entropy of a static 
black hole in MOG, due to thermal fluctuations. If we 
use the definition of the entropy from its temperature, 
we obtain 


/Sx — 7tG]\i Af^ (\/Of -|- 1 T 1 T cr) (\/or T 1 T 1) 

--In ^- - -- - - 

2 _(\/(a + l + l + (a)(\/(a + l + l) 


.( 78 ) 
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We also have 


A5 = 

= 7rGArM^(a(l + vT^f^) 

1 (1 + Of + \/\ + d) 

~2 (i + ^/TT^) 


(79) 


The corrections are the standard logarithmic correc¬ 
tions for the entropy of the MOG black hole. They re¬ 
duce to the entropy of a black hole expected from the 
holographic principle. Hence, these thermal fluctuations 
also lead to a violation of the holographic principle for 
these black holes. This is because the entropy of the black 
hole is reduced due to thermal fluctuations. This term 
by which the original entropy of the black hole is reduced 
is proportional to InS'oT^. It may be noted that as there 
are two different ways of defining MOG entropy i.e., Sto 
and Sao, the exact value of this term depends on such 
a definition chosen. However, in both these cases, the 
entropy of the MOG black hole reduces due to thermal 
fluctuations, and so these black holes will have less en¬ 
tropy than what would be expected from the holographic 
principle. Such violation of the holographic principle only 
occurs at small scale, where the temperature of the MOG 
black hole is sufficiently large, and so the effect of ther¬ 
mal fluctuations cannot be neglected. It may be noted 
that at this stage, the effects of quantum fluctuations 
can also not be neglected. In fact, the quantum fluctua¬ 
tions in the geometry of the black hole lead to these ther¬ 
mal fluctuations in the thermodynamics of the black hole. 
The violation of the holographic principle due to quan¬ 
tum fluctuations has been studied previously [35l|36]. We 
have demonstrated explicitly that such an effect occurs 
for MOG black holes due to thermal fluctuations. 


VII. CONCLUSIONS 

We have analyzed the thermodynamics of black holes 
in a modified theory of gravity that contains a repul¬ 
sive vector field with the gravitational charge Q = 
In the absence of matter, the field equations 
have a non-zero energy-momentum tensor formed from 
the field. For non-vanishing values of the param¬ 
eter a there are no pure vacuum solutions of the mod¬ 
ified gravitational field equations. The thermodynamic 
properties of the solutions corresponding to a spherically 
symmetric black hole and a rotating black hole were ana¬ 
lyzed, and it was found that the the entropy area law gets 
changed by increasing the size of the MOG parameter a. 
The thermodynamics of a regular black hole solution in 
MOG with a < acrit with two horizons was investigated. 
We also considered a regular solution without horizons 
describing a ”gray“ hole. There is a small but finite 
probability for particles to escape the gray holes. Thus, 
the modification of gravity which was initially proposed 
to explain the dynamics of galaxies and galaxy clusters 


without dark matter could end up predicting the absence 
of horizons for massive collapsed objects in our universe. 

We also analyzed the corrections to the thermodynam¬ 
ics of a black hole in MOG. The standard partition func¬ 
tion for a static black hole in the modified theory of grav¬ 
ity was used to perform this analysis. Thus, the thermo¬ 
dynamic properties of a MOG black hole were obtained 
using this partition function, and the thermal fluctua¬ 
tions to the thermodynamic quantities were computed. 
Motivated by a conformal field theory description, an ex¬ 
plicit form for the corrections to the entropy of the static 
black hole was found. It was demonstrated that correc¬ 
tions had the standard logarithmic form. The entropy 
of the black hole is less than what is expected from the 
holographic principle. Hence, the thermal fluctuations 
also lead to a violation of the holographic principle for 
static black holes in MOG. 

It has been observed that a traversable MOG wormhole 
exists as a solution of the MOG field equations [50]. So, 
it will be interesting to analyze the effects of the correc¬ 
tions to the MOG thermodynamics for such a wormhole 
solution. 

Corrections to ordinary GR black holes have been ob¬ 
tained using the generalized uncertainty principle m- 
It was demonstrated that such corrections can lead to 
the existence of black hole remnants, which can have im¬ 
portant phenomenological consequences Hg. It will be 
interesting to analyze such corrections for MOG black 
holes. Additional interesting and testable phenomeno¬ 
logical characteristics can be extracted from the MOG 
black hole solutions, including quasinormal modes and 
gravitational wave signatures of binary inspirals. The 
latter is of particular interest due to the definitive finger¬ 
prints of MOG in a black hole’s shadow m, which are 
expected to be detectable by the Event Horizon Telescope 
observations. If similar MOG effects are as pronounced 
for gravitational waves, then they should be detectable 
in the upcoming run of LIGO HU- Lastly, we can also 
consider other types of non-linear kinetic terms for MOG 
|52H58] 


JCi = {-BY, 

/C 2 = 462(1 - yi - (26)-iB), 

^3 = 6^(exp(-6“2B) - 1), 

/C 4 = -Sb^\n{l + b-^B), (80) 

where 5 , b are two parameters. It would be interesting 
to analyze the phenomenological consequences of such 
terms. 
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